Optical response of graphene under intense terahertz fields 
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Optical responses of graphene in the presence of intense circularly and linearly polarized terahertz 
fields are investigated based on the Floquet theory. We examine the energy spectrum and density 
of states. It is found that gaps open in the quasi-energy spectrum due to the single-photon/multi- 
photon resonances. These quasi-energy gaps are pronounced at small momentum, but decrease 
dramatically with the increase of momentum and finally tend to be closed when the momentum is 
large enough. Due to the contribution from the states at large momentum, the gaps in the density 
of states are effectively closed, in contrast to the prediction in the previous work by Oka and Aoki 
[Phys. Rev. B 79, 081406(R) (2009)]. We also investigate the optical conductivity for different field 
strengths and Fermi energies, and show the main features of the dynamical Franz-Keldysh effect in 
graphene. It is discovered that the optical conductivity exhibits a multi-step-like structure due to 
the sideband-modulated optical transition. It is also shown that dips appear at frequencies being the 
integer numbers of the applied terahertz field frequency in the case of low Fermi energy, originating 
from the quasi-energy gaps at small momentums. Moreover, under a circularly polarized terahertz 
field, we predict peaks in the middle of the "steps" and peaks induced by the contribution from the 
states around zero momentum in the optical conductivity. 

PACS numbers: 73.22.Pr, 78.67.Wj, 42.50.Hz 
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I. INTRODUCTION 

Since the experimental realization of graphene^ a 
monolayer of carbon atoms arranged in a honeycomb 
lattice, this material has aroused enormous interest due 
to its unique physical characteristics^— Among differ- 
ent works in this field, the linear optical property of 
graphene is one of the main focuses of attention*^— The 
theoretical works based on the Dirac Hamiltonian show 
that the optical conductivity at high frequency is domi- 
nated by the interband optical conductivity, which takes 
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a constant value of ao = for frequency larger than 
twice of the Fermi energy Ep and approaches to zero 
for frequency below 2Ef due to the Pauli blocking.— ~— 
The universal value <tq of optical conductivity has been 
demonstrated to be valid not only in the noninteracting 
limit , 12 ' 13 but also in the presence of the disorder and the 
electron-electron interaction as long as the Dirac-cone 
approximation remains validJ^— The constant optical 
conductivity in a wide range of frequency also has been 
observed in the optical experiments*^— On the other 
hand, the optical conductivity at low frequency is de- 
termined by the intraband optical conductivity, which 
presents a Drude peak centred at zero frequency and is 
strongly influenced by the sample- dependent scattering 
behavior J^r— 

Recently, influence of an intense ac field on electri- 
cal and optical properties in graphene has also attracted 
much attention.— ~— It has been found that the ap- 
plication of an intense ac field can dramatically mod- 
ify the band structure and hence the density of states 
(DOS).— ~— Their results also showed that a stationary 
energy gap appears around the Dirac point under a cir- 



cularly polarized ac fiekh 21 ' 24 Oka and Aok i 21 ' 22 i 26 calcu- 
lated the dc and ac conductivities in graphene irradiated 
by an intense circularly polarized light via the extended 
Kubo formula, and proposed the photovoltaic Hall ef- 
fect, which is a novel Hall effect occuring in the absence 
of uniform magnetic fields. The dc transport properties 
of graphene-based p-n junctions under an intense ac field 
were also investigated theoretically! 20 ' 27 However, contri- 
bution in the optical spectra in graphene from the optical 
sidebands^ 7 . - — has not yet been well investigated. 

In semiconductors, the contribution from the side- 
bands has been demonstrated to be important for the 
optical and transport properties. Many interesting 
phenomena, such as the photon-assisted tunneling , 48 ' 49 
the sideband generation of exciton^ - — the ac Stark 
effect^i 4 !^ and the dynamical Franz-Keldysh (DFK) 
effect p 3 - - — as well as spin generation and manipulation 
utilized by the intense terahertz (THz) fielder— are re- 
lated to the formation of sidebands. Among these ef- 
fects, the DFK effect describes the influence on optical 
spectra from the sidebands of the expanded states , 43 ' 44 
which includes finite absorption below the band edge 
from the contribution of the sidebands below the bottom 
of the conduction band and the step-like behavior above 
the band edge due to the sideband-modulated general- 
ized DOS. Just as the DFK effect in semiconductors, the 
formation of optical sidebands should also influence the 
optical spectra near the absorption edge around 2Ep in 
graphene. Nevertheless, the band structure of graphene 
is gapless and the energy dispersion is linear, which is 
quite distinct from semiconductors. Thus the DFK ef- 
fect in graphene is expected to present some unique be- 
haviors. This makes the investigation on this problem 
become very interesting and important. It is also noted 
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that in previous investigations 6 , only the optical conduc- 
tivity with Ep much smaller than the frequency of the 
ac field was discussed and thus the contribution from the 
optical sidebands is difficult to identify. In the present 
work, we calculate the optical conductivity of graphene 
under the intense THz field for various Fermi energies 
and field strengths in order to gain a complete view of 
the DFK effect in graphene. 

In order to include the contribution from the opti- 
cal sidebands explicitly, we solve the time-dependent 
Schrodinger equation by using the Floquct theory 55 
and obtain the optical conductivity in graphene un- 
der an intense THz field via the nonequlibrium Green 
functions ! 44 ' 56 In this paper, we focus on the optical 
conductivity at high frequency, which is known to be 
insensitive to the scattering strength»i£r— This allows 
us to ignore the detail of the scattering and only dis- 
cuss the optical conductivity in the noninteracting limit. 
We first examine the energy spectrum and DOS. It is 
found that gaps appear in the quasi-energy spectrum due 
to the single-photon/multi-photon resonances ! 41 ' 48 These 
quasi-energy gaps are pronounced at small momentum, in 
consistence with the previous investigations! 20 ' 21 ^ 4 How- 
ever, the quasi-energy gaps decrease dramatically with 
the increase of momentum and finally disappear when 
the momentum is large enough. Therefore after taking 
account of the contribution from the states with large 
momentum, gaps in the DOS are effectively closed, in 
contrast to the prediction by Oka and Aokii^i Our re- 
sults of the optical conductivity reveal the main features 
of the DFK effect in graphene. 

This paper is organized as follows. In Sec. IIA, we ob- 
tain the energy spectrum by exploiting the Floquct the- 
ory. Then in Sec. IIB, we derive the optical conductivity 
via the nonequlibrium Green functions. The numerical 
results of the DOS and optical conductivity are presented 
in Sec. III. Finally, we summarize in Sec. IV. 

II. MODEL AND FORMALISM 

A. Hamiltonian 

We consider a graphene layer placed in the x-y plane. 
In the vicinity of the Dirac point, the effective Hamilto- 
nian of graphene can be written as (ft = I) 5 -^ 

i?o(k) = v F (fia x k x + a y k y ). (I) 

Here // = 1(— 1) for K(K') valley; vp is the Fermi veloc- 
ity; k represents the two-dimensional wave vector rela- 
tive to K(K') point; a is the Pauli matrix in the pseu- 
dospin space formed by the A and B sublattices of the 
honeycomb lattice. Here and hereafter, symbols with " 
present the 2x2 matrices in the pseudospin space. The 
eigenvalue and eigenvector of Hq are E^ v — vv-p\k\ and 
£ ki/ = 1/V2(five ^ l t l6k 1 1) T , respectively, with v being 1 
(— I) for electron (hole) band and # k representing the po- 
lar angle of k. Substituting k by k + eA(t), one obtains 



the effective Hamiltonian in the presence of a THz field 

£&(k,t) - HSM+H^it), (2) 
#thz(*) = ev F [fia x A x (t) + a y A y (t)}. (3) 

For convenience, we choose the THz field as E(t) = 
Eo(e x cos#e cos fit + e y sin#E sin fit). Thus the vec- 
tor potential A(t) reads A(t) = e x cos6*e sinfii + 
e y sin0 E cos fit). For E = 0, tt/4 and tt/2, the THz fields 
are linearly polarized along the x-axis, circularly polar- 
ized and linearly polarized along the j/-axis, respectively. 
Without loss of generality, we set the THz field linearly 
polarized along the x-axis when discussing the case for a 
linearly polarized THz field. 

By exploiting the Floquet theory; 5 - 5 - the solution of the 
Schrodinger equation id t \$£ a (t)) = fl^(k, t)\$£ a (t)) 
has the form ( |3> ka (i)) 1S referred to as the Floquet state 
in the following): 

oo 

\K a (t)) = e- i£ ^ t £ e^|^ kQ ), (4) 

n— — oo 

in which a represents the branch index of the solu- 
tion; ££ Q and |V>" kQ ) = (^ka.^k«) T are the eigenvalue 
(quasi-energy) and eigenvector determined by 

" = y[t>cos0 E - <r S ine E )^-" 

- (/x cos $e +^^^"1 + ^Me-^ ek < k - CT -(5) 

with (3 = v-peEo/fl' 2 . The above equation shows that the 
relation between the normalized quasi-energy and 
the normalized momentum upk/fi is only determined by 
the dimensionless quantity j3£& 

Due to the periodicity of H^ s (k, t), the eigenvalues are 
also periodic, i.e., if Et is a solution of Eq. ©, then 

^ka = ^lia + ^ i s also a solution. It is evident that the 
eigenvectors of and ££ a satisfy ip 1 ^ = ip^L" > tllus 
^k'a an d ^ka correspond to the same physical solution 
of the Schrodinger equation. Namely, the quasi-energy 
is a multi-valued quantity of the Floquet stated Never- 
theless, for each momentum and each valley, the number 
of the independent quasi-energies is 2, which is deter- 
mined by the dimension of the Hilbcrt space. For conve- 
nience, we choose the independent quasi-energies e k?? in 
the reduced Floquet zone (— f2/2, f2/2], with ij = ± repre- 
senting the index of independent solutions. These quasi- 
energies are referred to as the reduced quasi-energies in 
the following. The corresponding eigenvectors are la- 
belled as |</>™ k)7 ). Therefore, by choosing the proper inte- 
ger I, arbitrary quasi-energy £ ka and the corresponding 
eigenvectors ip^ a can be written into the form 

£L=e£n + M, (6) 
& = ^4 CT - (7) 
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In addition, for the reduced quasi-energies with different 
r), one has 



-k+ 



y l Tjuk- 



(8) 
(9) 



From Eq. (j4j), it is seen that the Floquet state |$k„(t)) 
(the general solution index a has been replaced by the 
independent solution index rj) contains Fourier compo- 
nents with different frequencies, quite distinct from the 
eigenstate of Hq, which only takes one eigen- frequency. 
Each Fourier component corresponds to a sideband of 
this Floquet state. For the Floquet states with reduced 
quasi-energy ej* , sidebands appear at the frequencies 
(quasi-energies) 



?k»7 fc k»7 

and the corresponding weights are 



(10) 



(11) 



From Eqs. © and ([§]), one can see that the quasi-energies 
and weights of the sidebands satisfy 



c 1 = w£r n . 



(12) 
(13) 



Besides the quasi-energy, another important quantity 
of the Floquet state is the mean energ y 28 i 59 ~ — 



1 



Tier* JJW7 
[Mk.7] t filer) 1 



(14) 



where Tq = 2n/£l is the period of the applied THz field. 
Independent of the choice of the quasi-energy, the mean 
energy is a single-valued structure quantity of the Flo- 
quet state. Thus in previous workS) 28 i 59 ' 60 the mean en- 
ergy was utilized to identify the filled Floquet states, i.e. 
states with lower mean energy will be occupied at first. 
In this paper, we also use this ansatz to obtain the dis- 
tribution function of the Floquet state, which will be 
discussed in the next subsection. Moreover, the mean 
energy is used to identify whether the Floquet state is 
electron- or hole-like. Analogous to the definition of the 
electron and hole states in the field-free case, we define 
the quasi-electron (ry = +) and quasi-hole states (r) = — ) 
as the Floquet states satisfying e£ + > and < 0, 
respectively. 



B. Optical conductivity 

It is known that the optical absorption is measured 
by the real part of optical conductivity*^— Therefore 
we focus on the real part of the optical conductivity in 



the following. For the probing light field of frequency 
toi with the polarization in the I (= x, y) direction, the 
linear-response theory yields the real part of optical con- 
ductivity: 



Rea u {T, uji) 



9s9v 



IniH^T,*,). 



(15) 



Here g v = 2 and g s = 2 are the valley and spin degenera- 
cies, respectively; n^(T, uj{) is retarded current-current 
correlation function in the K valley. Here and hereafter 
we only give the correlation function in the K valley and 
omit the valley index /i in all symbols, as the contribu- 
tions to optical conductivity from both valleys are iden- 
tical. n[ ; (T, uji) can be written as 

/oo 
Jre^WuiT+^T--), (16) 



where 



W u {t,t') = -i6{t-t'){\j l {t),f{t')]) (17) 



with y = evpai presenting the I component of the current 
operator in graphene . 12 ^ 14 Via the nonequlibrium Green 
function method i 44 i 56 we have 



nFz(t,f) 



:cr 1 cr2 



+ 3L in ^w; (*> t')jL^ 2 (f, t)] , (18) 



with G^'^\t,t') representing the retarded (advanced, 
lesser) single-particle Green function^ Substituting 
Eq. (32]) into Eq. lfl6"]>. one has 



nS(T, Wl ) 



E 



2^ 



'jka 1 a 2 jktj' 1 a' 2 G^ cr ' 2<7l (T, U)) 



x [GL 2CT ; (T, lo + WI ) + Gl a2< (T, u - w ,)] .(19) 

S^ce [G<^ (t,t>)]* = -G^^t) and 

[GL 1CT2 (M')]* - Gi a2ai (t',t), W u {t,t') is real. Thus 
one obtains 



ImH[ ; (T,^) = -\ir u (T, Ul ) - n[ ; (T, — w/)] 



E 



do; 
2^ 



A kCT2CT; (T, w -^)]G< ffi (7». 



(20) 



In above equations, we have used the relation 
A kffl(T2 (7» = i[Gi aia2 (T,u;) - G^ ia2 (T,u,)}. In the 
noninteracting limit, the retarded Green function is given 
by 

Gl(t 1 ,t 2 ) = ~ t 6(t 1 -t 2 )J2\$k V (t 1 ))(^(t 2 )\. (21) 
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Thus the spectral function in the frequency space can be 
written as 

/oo 
dre^£|$ k „(r+l)><$ k „(T-I)| 
-°° n z z 

= 2tt ]T e i( "-" l)OT |0k„>(C^(^ - £k, + (n 

nmr} 

+ m)|). (22) 

The next step is to calculate (T, w). The equal-time 
lesser function can be expressed in the form 

=iJ2t*™M^vM(^nM- (23) 

Here p kr?ir)2 (t) is the nonequilibrium density matrix in the 
Floquet picture ) 54 ! 62 which can be determined by the ki- 
netic equation including the electron-impurity, electron- 
phonon and electron-electron scatterings * 54 ' 56 ! 62 ' 63 This 
approach is very complicated and left as the subject of 
our future work. Here we obtain the lesser function based 
on a simple ansatz following the previous works; 28 ' 59 ' 60 
the steady-state density matrix is diagonal with the diag- 
onal term / k?7 being the Fermi distribution on the mean 
energy of the Floquet state, 

Pkrn^it) = far)5r nn2 = n F (^kri)S m rj 2 • (24) 

In the present paper, we focus on the case at zero tem- 
perature, and then 

h v = 0{e kv -E F ), (25) 
where the Fermi energy E-p can be obtained from 

Ne = 9s9v ( E V - £kr,)- (26) 

Substituting Eq. ((24]) into Eq. ([23]>. one gets 

G^(t,t)=iYf^Ur,(t))(^r,(t)\- (27) 
n 

By further exploiting the generalized Kadanoff-Baym 
Ansatz^ 6 - the two-time lesser function can be obtained 
as 

G<(t,t') = iGi(t,t')G<(t',t')-iG<(t,t)GUt,t') 

= ^/k^ kl7 (t)X<M*')l- (28) 
n 

Then we have 

G<(T,u l )=iJ2hr,A kf ,(T,uj), (29) 
n 



with i k „eZ» = J^dre— |$ k „(T+f))($ k ^T-§)|. 
Using the above equations, one obtains the optical con- 
ductivity 

Rea u (T,^) = ^ £ (C^ilC^d \&\<V 

m 1 Ti 2 m 2 

ft 

- £k V2 - (»i +mi - n 2 - Tti2 ) ~ + wi).(30) 
The time-averaged optical conductivity can be written as 

Reafpiuii) = ±- [ ° dT Reoi,(T,w,) 
J o Jo 

x (A m - /kr )2 )<5(e k ,, 1 - e kl)2 - Nfl + w{). (31) 

It is found that the time-averaged optical conductivity 
from the nonequlibrium Green functions agrees with that 
from the extended Kubo formula used in Ref. [26| (see 
also Appendix |A"| . However, it is noted that our method 
not only gives the time-averaged optical conductivity but 
also the time-dependent one, and thus can provide the 
dynamical information of optical response, which can 
be observed via the time-resolved measurements. It is 
also noted the distribution used by Oka and Aoki,— i.e., 
/k£ = EM-EWjKCki^L,)! 2 ' is obtained by projecting 

the equilibrium Fermi distribution in the field-free case 
to the basis set formed by the Floquet eigenvectors, cor- 
responding to the "sudden approximation". Obviously, 
this distribution, referred to as the projected distribution 
in the following, is quite different from the mean-energy- 
determined distribution used in this paper. The compar- 
ison of the optical conductivities obtained from these two 
distributions will be addressed in next section. 



III. NUMERICAL RESULTS 

In this section, we discuss the numerical results of 
the energy spectrum, DOS and optical conductivity of 
graphene under an intense THz field. The typical pa- 
rameters used in the computation are ft = 5 THz, 
E = 15 kV/cm (corresponding to j3 — 2.3), N e = 1.5 x 
10 12 cm" 2 (E F = 60) and 2.5 x 10 11 cm" 2 (E F = 2.7ft). 
It is noted that our results can be generalized to the other 
ft regime as long as Eq. (|2~4")l remains valid, since the be- 
haviors of the energy spectrum and the DOS are only 
determined by (3 and the behavior of the optical con- 
ductivity is only determined by /3 and E F /fl£^ We also 
restrict our investigation of optical conductivity in the 
frequency regime toi > ft, since the scattering process is 
not considered in our model and the optical conductivity 
at low frequency is known to be strongly dependent on 
the scattering strength^ - — 
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FIG. 1: (Color online) Circularly polarized THz field with 
f3 — 0.4. (a) Quasi-energies f{L against the normalized mo- 
mentum. The color coding represents the weight W\L of the 
corresponding sideband. The blue crosses represent the se- 
lected quasi-energies of the quasi-electron states, (b) Quasi- 
energy £k+ (red solid curve) and mean energy £k+ (blue 
dashed curve) of the quasi-electron state as well as the field- 
free electron energy (green dotted curve) against the normal- 
ized momentum. In the inset of (b), we plot the quasi-energy 
around the gap at k = 0.5f2/«F and e-^ v = 0.557 [the re- 
gion labelled by the box in (a)] as well as the corresponding 
mean energy as function of momentum. The solid (dashed) 
curves with and without crosses represent the quasi-energies 
(mean energies) of the quasi-electron (r/ = +) and quasi-hole 
(77 = — ) states, respectively. The black chain line indicates 
the momentum of the crossover point ko- Note the scale of 
mean energy is on the right-hand side of the frame. 



A. Energy spectrum 

Although the energy spectrum in this system has been 
investigated by many works^2r— a complete investiga- 
tion on this problem is still lacking. In particular, the 
energy spectrum at large momentum has not been well 
investigated in previous works. In this section, we inves- 
tigate the energy spectrum under a circularly polarized 
THz field in a wide range of momentum. Pronounced 
quasi-energy gaps appear at small momentum in both 
cases with low and high field strengths. We also discuss 
the case with a linearly polarized THz field and show that 
the energy spectrum becomes anisotropic. 



FIG. 2: (Color online) Circularly polarized THz field with 
/? = 2.3. (a) Quasi-energies fjL against the normalized mo- 
mentum. The color coding represents the weight W{L of the 
corresponding sideband. The blue crosses represent the se- 
lected quasi-energies of the quasi-electron states, (b) Quasi- 
energy Ek+ (red solid curve) and mean energy Ek+ (blue 
dashed curve) of the quasi-electron state as well as the field- 
free electron energy (green dotted curve) versus the normal- 
ized momentum. The thin black lines mark the mean ener- 
gies being 1.557, 1.857, 2.357 and 2.757, which are used in the 
Fig. [6jc). Note the scale of mean energy is on the right-hand 
side of the frame. 



We first concentrate on the case for a circularly po- 
larized THz field with low field strength (small /3). 
In Fig. QJa), the quasi-energies of the sidebands £]L 
[Eq. (fTU)) ] for j3 — 0.4 are plotted as function of mo- 
mentum. Here the color coding represents the weight 
W£ [Eq. (fTTj) ] of the corresponding sideband. Since 
the quasi-energy spectrum is isotropic under a circu- 
larly polarized field, we only present the results with 
6>k = 0. The most interesting feature seen in this fig- 
ure is the appearance of gaps at small momentum in 
the quasi-energy spectrum, in consistence with the previ- 
ous investigations! 20 ' 21 ' 24 These gaps appear around the 
momentums toS!/2«f and the quasi-energies m£l/2 + IQ, 
with m and I being integers. All gaps at the same mo- 
mentum share the identical magnitude due to the period- 
icity of the quasi-energy spectrum. These quasi-energy 
gaps can be attributed to the ac Stark splittings induced 
by the single-photon/multi-photon resonances! 41 ' 48 ' 59- — 
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FIG. 3: (Color online) Linearly polarized THz field with f3 = 2.3. Quasi-energies against the normalized momentum for 
different polar angles #k- The color coding represents the weight W\L of the corresponding sideband. fc„ in (c) stands for k ■ n 
with n being the unit vector at an angle 7r/6 with respect to the x-axis. 



The physics is that if a pair of states are coupled by an 
electromagnetic field with frequency f2 and the energy 
difference between these two states equals to mil, then 
an ac Stark splitting appears unless the corresponding 
transition is forbidden. Specifically, in the present case, 
the gaps around the momentum satisfying 2i>pfc = mil 
are induced by the m-photon resonances^ 

Nevertheless, the behaviour of the quasi-energy spec- 
trum in large momentum regime becomes very different. 
From Fig. QJa), one can see that the energy gaps de- 
crease dramatically with the increase of k and finally 
tend to be closed when the momentum is large enough. 
This effect can be understood via Eq. (IB24[) under the 
rotating-wavc approximation (the derivation is presented 
in Appendix [Bj, which shows that the quasi-energy gap 
around the momentum mfl/2vF is determined by the ef- 
fective coupling parameter |y m (#k)|; given by [see also 
Eq. (S3)] 

MftOI = -j=\J m +i(y/20) - J TO _i(\/2/3)|, (32) 

with J m (y/2f3) being the Bessel function. Due to the in- 
dependence of J m (a/2/3), the effective coupling parame- 
ter |y m (0k)| decreases dramatically with the increase of 
m when m is large enough. Thus the quasi-energy gaps 
become negligible at large momentum. Our calculations 
also show that, associated with the absence of the gaps, 
the band- mixing, i.e., the hole (electron) component of 
the quasi-electron (quasi-hole) state, becomes negligible, 
and the wave function becomes very close to the one with- 
out the interband term of Htuz ■ 

In Fig. Hlb), we plot the mean energy £k+ [Eq. (|T4| ] 
and the quasi-energy £k+ of the quasi-electron state 
(£k+ > 0) as well as the field-free electron energy for 
/3 = 0.4 as function of momentum. In order to compare 
the above three energies directly, we choose the quasi- 
energy £k+ according to the following rules: 1) the conti- 



nuity of the quasi-energy is kept as far as possible; 2) the 
quasi-energy at large momentum is closest to the field- 
free electron energy. We also plot these selected quasi- 
energy in Fig. [TJa) (blue crosses). From Fig. [TJb), it is 
shown that the mean energy and the quasi-energy are far 
away from (very close to) the field-free energy at small 
(large) momentum, in consistence with the behaviour of 
the quasi-energy gap. Moreover, it is also seen that the 
mean energy reaches zero at a momentum somewhere 
inside the quasi-energy gap, except the one at zero mo- 
mentum. In order to reveal the underlying physics, we 
plot the quasi-energy around the gap at k = 0.5Q/vf and 
fkr; = 0.5f2 [the region labelled by the box in Fig.QTa)] as 
well as the corresponding mean energy as function of mo- 
mentum in the inset of Fig. [IJb). It shows that there is 
a crossover point (labelled as k ) between the quasi-hole 
state (blue solid curve) and the quasi-electron state (red 
solid curve with crosses) in the quasi-energy above the 
gap. Since the quasi-energy varies continuously with k 
at this point, the corresponding mean energy also varies 
continuously, i.e., £k + = £k -- Also, from Eqs. ([8]) and 
©, one can see that £k + = — f?k -- Consequently, the 
mean energies at this point must be zero. The finite mean 
energy at k = is based on the similar reason. As shown 
in Fig. []Ja), the quasi-energy above the gap at zero mo- 
mentum always belongs to the quasi-electron state [blue 
crosses in Fig. HJa)], hence no crossover appears. 

Now we turn to the case for a circularly polarized THz 
field with high field strength (Fig. [2]). From Fig. 0a), 
it is seen that pronounced quasi-energy gaps appear 
in a wider range of k. However, unlike the previous 
case, the momentums of the gaps markedly deviate from 
toO/(2wf)- This is because the effect of the rapidly 
varying terms in the resonance equations [Eqs. (IB15P 
and (|B16[) ] cannot be neglected for strong field, i.e., the 
rotating-wavc approximation is not valid. The joint effect 
of the terms with different oscillating frequencies leads to 
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the complicated behaviour of the quasi-energy spectrum. 
However, the quasi-energies where the gaps appear are 
still around ifi/2, determined by the symmetry of the 
quasi-energies of the sidebands [Eq. (|T5|) ]. Figure [5] also 
shows that the quasi-energy gaps become extremely small 
and the quasi-energy and mean energy become very close 
to the field-free energy at large momentum. The physics 
is similar to the weak-field case. When the photon num- 
ber involved in the resonance becomes too large, the ef- 
fective resonant coupling becomes extremely weak. Thus 
the influence of the THz field on the energy spectrum 
becomes negligible. 

Finally we address the case for a linearly polarized THz 
field (Fig. [3]). Here we only present the results with high 
field strength as they are sufficient for showing the main 
features in this case. Recall that the linearly polarized 
THz field is set along the cc-axis throughout this paper, 
thus #k equals to the angle between the momentum and 
the THz field. From Fig. [3j it is seen that quasi-energy 
gaps are closed at zero momentum. This effect can be 
understood from the exact analytical solution 



|<&k=o±(*)) = £j±»G8)e inn *(l, ±1) 



(33) 



Clearly, the quasi-energy is zero at zero momentum, thus 
the gap disappears. Another interesting feature is the 
anisotropy of the quasi-energy spectrum^ For momen- 
tum along the x-axis [Fig.[3ja)], all gaps are closed, since 
the interband term of Hthz [Eq. (|B4[) ] becomes zero 
and hence the quasi-energy is exactly the same as the 
field-free energy, as shown in Eqs. (IB5|) and (|B6[) in Ap- 
pendix [HI For momentum along the y-axis [Fig. EJb)], 
all gaps except the ones around k — 0.5il/vp are ef- 
fectively closed. This is because the intraband term of 
Hthz [Eq. (|B4I) ] becomes zero, and thus the gaps from 
the multi-photon resonances become extremely small, as 
shown in Fig. [5] in Appendix [B] For other polar angle, 
e.g., #k = tt/6 [Fig. [3jc)], pronounced quasi-energy gaps 
appear not only around k = O.bfl/v-p but also around the 
other resonant points at small momentum, similar to the 
case with a circularly polarized THz field. 



B. DOS 

Then we turn to investigate the DOS. By using the 
spectral function Eq. (f2"2"j) . one obtains the time- averaged 
DOS 



D avc {uj) 



9s9v 

27I"To Jq 



dT ^Tr4,(T,w) 



9s9, 



£«JV£aW W - £ka)- (34) 



ko 



The last equality is derived from Eqs. ^ and J?]). It 
is noted that the above equation is in the same form 
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FIG. 4: (Color online) (a) Time-averaged DOS versus u) 
under circularly (red solid curve) and linearly (blue dashed 
curve) polarized THz fields for /3 = 2.3. The region close to 
oj ~ is enlarged in the inset, (b) Time-averaged DOS ver- 
sus lo under a circularly polarized THz field for ft = 2.3 with 
all relevant states (red solid curve), limited in the momen- 
tum regimes v?k < 2Q. (blue dashed curve) and vf& < 0.5fi 
(green dotted curve). The yellow chain curve represents the 
DOS from Ref. EH. 



as the one reported by Oka and Aoki.— In Fig. |4ja), we 
plot the time-averaged DOS under circularly and linearly 
polarized THz fields for (3 = 2.3. (Note that j3 = 2.3 cor- 
responds to F/n = = 1.6 in Ref.HH^) The region 
close to to ~ is enlarged in the inset. Here we only 
show the DOS in the positive frequency regime (elec- 
tron regime), since the DOS in the negative regime (hole 
regime) is symmetrical to the positive one. 

We first focus on the case for a circularly polarized THz 
field [red solid curve in Fig. Hta)]. It is seen that the gaps 
in the DOS are effectively closed, in contrast to the pre- 
vious report by Oka and Aoki<2i The underlying physics 
is as follows: as shown in Fig. [21 a), the quasi-energy gaps 
are closed when the momentum is large enough. The con- 
tribution from these states closes the gaps in the DOS. 
We also plot the DOS with all related states (red solid 
curve), limited in the momentum regimes vpk < 2f2 (blue 
dashed curve) and v^k < 0.50 (green dotted curve) in 
Fig.HKb). :t is shown that the DOS with v F k < 2ft (blue 
dashed curve) is almost the same as the one shown in 
Ref. [2l] (yellow chain curve) . This indicates that the con- 
tribution from large momentum is not properly counted 
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FIG. 5: (Color online) Circularly polarized THz field, (a) 
Time-dependent optical conductivity as function of the op- 
tical frequency and time with /3 — 2.3. (b) Time-averaged 
optical conductivity as function of the optical frequency with 
different field strengths. Ep = 6fi in both figures. 

in the previous investigation^ 1 In addition, it is seen 
that sharp peaks appear in the DOS. This effect origi- 
nates from the van Hoff singularities which appear at the 
momentum of all quasi-energy gaps due to the isotropic 
quasi-energy spectrum^ 

Then we turn to the case of a linearly polarized THz 
field [blue dashed curve in Fig.[4|a)]. It is also shown that 
the gaps are absent. However the physics is different here. 
As shown in Fig. [S^a), the quasi-energy gaps disappear 
for the momentum along the a;-axis, even at small k, and 
the gaps due to the multi-photon resonances become neg- 
ligible for momentum along the y-axis. Therefore there 
is no gap in the DOS. Moreover, the peaks from the van 
Hoff singularities are less pronounced due to the different 
nature of the van Hoff singularities here associated with 
anisotropic quasi-energy spectrum.— 



C. Optical conductivity 

1. Under a circularly polarized THz field 

In this subsection we discuss the optical conductivity 
under a circularly polarized THz field. Without loss of 



1,4 





1.2 




i 


E F = 2.7 n 

2.3 n 








■ ■ 


1.8 £2 





1 






i.5 n 


> g 


0.8 








V 










o 

a. 


0.6 








0.4 
0.2 


-' 


1 'r v, 


'If' f if 



1 

12345 6 789 

(0,/n 



FIG. 6: (Color online) Circularly polarized THz field, (a) 
Time-dependent optical conductivity as function of the op- 
tical frequency and time with /3 = 2.3 and Ef = 2.7Q,. (b) 
Time-averaged optical conductivity as function of the opti- 
cal frequency with different field strengths. Ep = 2.7J7. (c) 
Time-averaged optical conductivity as function of the optical 
frequency with /3 = 2.3 for different Fermi energies. 



generality, we restrict ourselves to the rt-type case, i.e., 
Ep > 0. Thus only the quasi-electron states with eh v > 
E-p are not occupied and the the corresponding interband 
transitions are allowed. 

We first focus on the case with high Fermi energy. The 
time-dependent and time-averaged optical conductivities 
as function of the optical frequency loi are plotted in 
Figs.jnja) and (b), respectively. Here Co = fw is used as 
the unit of the optical conductivity. It is shown that the 
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optical conductivity exhibits a multi-step-like structure 
at u)i ~ 2Ep, in contrast to the single-step-like behaviour 
in the field- free case [yellow chain curve in Fig. EJb)]. 10-19 
Similar to the DFK effect in semiconductors^ this ef- 
fect is from the sideband-modulated optical transition, 
i.e., nfi in the delta functions in Eqs. (|30|) and ([31]) . It 
is also seen that the number of "step" increases with the 
increase of the field strength. This can be understood by 
noticing that the weight of the sideband is distributed in 
a wider range of frequency for stronger field. 

Figure [5] also shows that the optical conductivity varies 
mildly with the increase of u>i in each "step" , which is 
quite different from the DFK effect in semiconductors 44 
where the optical absorption is strongly dependent on oji 
in each "step" . This behavior can be understood as fol- 
lows. Since the effect from the interband term of -Hthz 
becomes negligible at large momentum, the optical con- 
ductivity at high Fermi energy can be approximately de- 
scribed by Eq. (|C5[) (the derivation is presented in Ap- 
pendix [C]), which indicates that the frequency depen- 
dence of the optical conductivity is only from the factor 
1 — Ntt/uji in the optical transition between the sidebands 
with the energy difference 2?k+ + Nfl. This factor orig- 
inates from the linear dispersion of graphene. It is also 
noted that the optical conductivity is dominated by the 
optical transitions with small \N\, and the pronounced 
"steps" only appear at u>i ~ 2Ep. Thus the frequency de- 
pendence of the optical conductivity becomes very weak 
in each "step" for high Fermi energy. 

Then we turn to the case with low Fermi energy. The 
time-dependent and time-averaged optical conductivities 
are plotted as function of the optical frequency in the 
cases with Ep — 2.7H in Figs. [UJa) and (b), respec- 
tively. It is seen that dips appear around the frequen- 
cies satisfying Ifl when the applied THz field is strong 
enough. This is because the states at small momen- 
tum, where the quasi-energy gaps become pronounced as 
shown in Fig. HJb), can contribute to the optical conduc- 
tivity in this case. More interesting features are presented 
in the Fermi energy dependence of optical conductivity 
[Fig-Hl c )]- It is shown that peaks appear in the middle of 
the "steps" for Ep = 2.3fl (red solid curve). The scenario 
is as follows. From Fig. [2Kb), it is seen that the maximum 
of the mean energy £k+ at 2.5fl/vp is slightly higher than 
the Fermi energy 2.30. Consequently, in this momentum 
regime only the quasi-electron states around 2.5Q/vp can 
contribute to the optical conductivity and hence induces 
the peaks in the middle of the "steps" . When the Fermi 
energy decreases, more and more states in this momen- 
tum regime can contribute to the optical conductivity. 
Consequently the peaks become wider and wider, and 
finally become the new "steps" , as shown in the case 
with Ep — 1.8S1 (yellow chain curve). Moreover, when 
the Fermi energy decreases to 1.5fl (blue dashed curve), 
sharp peaks appear in the optical conductivity. This orig- 
inates from the contribution around k = 0, in which £]£+ 
is slightly higher than 1.5f2 as shown in Fig. HJb). One 
also notes that the van Hoff singularities at nonzero mo- 
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FIG. 7: (Color online) Linearly polarized THz field with 
/3 — 2.3. (a) Time-averaged optical conductivity as func- 
tion of the optical frequency for different Fermi energies, 
(b) Time-averaged optical conductivities versus the optical 
frequency from the states in different polar angle regions 
[0k - 7r/128,0k + tt/128) as well as the one with all polar 
angles integrated. Ep = 2.5Q. (c) Mean energies of the 
quasi-electron states with different polar angles 6^ against 
the normalized momentum. The thin black line indicates the 
mean energy being 2.5f2. 



mcntum have no effect on the optical conductivity with 
finite Fermi energy, since the mean energy is close to zero 
at the momentum of the quasi-energy gaps, as shown in 
Fig. (Kb). 
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2. Under a linearly polarized THz field 

We also plot the time-averaged optical conductivity as 
function of the optical frequency under a linearly polar- 
ized THz held in Fig. EJa). The behavior in this case is 
simpler than that under a circularly polarized THz field. 
The multi-step-like behavior and the dips around iO still 
appear. Nevertheless, the peaks in the middle of the 
"steps" , which still appear in the optical conductivity at 
low frequency for E-p = 2.50 and 1.2f2, are much less pro- 
nounced than the ones under a circularly polarized THz 
held. In order to reveal the underlying physics, we plot 
the time-averaged optical conductivity from the states 
with different polar angles as function of the optical fre- 
quency for E-p — 2. 557 in Fig. [7Jb). The corresponding 
mean energies are plotted against the normalized mo- 
mentum in Fig. [7]Jc). As mentioned above, the peaks in 
the middle of the "steps" only appear in the situation 
satisfying the criteria that a local maximum of the mean 
energy £k+ is slightly higher than the Fermi energy. From 
Fig. Etc), one can see that the above criteria is satished 
for the states with 9^ = ir/12 and 7r/6, thus pronounced 
peaks appear in the optical conductivity from the states 
with these angles [red solid and blue dashed curves in 
Fig. [7Jb)] and the corresponding angles with the same 
quasi-energy^ However, this criteria cannot be satished 
for the states with the other polar angles, e.g., 8^ = n/8 
and 0k = 7r/4, so peaks are absent in the correspond- 
ing optical conductivity [green dotted and yellow chain 
curves in Fig. E{b)]. Therefore, after the summation of 
the contribution from all polar angles, the peaks in the 
middle of the "steps" become much less pronounced. 

3. Comparison of the optical conductivities calculated with 
the mean-energy-determined distribution and the projected 
distribution 

In this subsection, we compare the optical conductiv- 
ities calculated with the distribution determined by the 
mean energy and the projected distribution. As men- 
tioned above, the projected distribution used in Ref. 
is described by (see also Appendix [AJ 

f^ = J2 n ^ E ^^v)\ 2 - ( 35 ) 

Recall that E^ and Ck^ are the eigenvalue and eigen- 
vector of Ho. From Eq. (pJ5j) . one can recognize the 
main features of the projected distribution at zero tem- 
perature: at k < kp = Ep/vp, both the field-free elec- 
tron and hole states are occupied, thus the corresponding 
quasi-electron and quasi-hole states are both occupied; at 
k > kp, only the hole state is occupied, so the distribution 
of the Floquet state is determined by the hole component 
of this state. Therefore one can see that only the states 
with k > kp can contribute to the optical conductivity, 
and the contribution decreases with the increase of the 
band-mixing. 
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FIG. 8: (Color online) Time-averaged optical conductivities 
from the mean-energy-determined distribution (solid curves) 
and the projected distribution (dashed curves) are plotted as 
function of the optical frequency under a circularly polarized 
THz field with f3 = 2.3 for different Fermi energies. 

The time-averaged optical conductivities from the 
mean-energy-determined distribution and the projected 
distribution are plotted as function of the optical fre- 
quency in Fig. |SJ It is seen that in the case with high 
Fermi energy, the results from both distributions are al- 
most the same. This is because the optical conductivity 
in this situation is determined by the contribution from 
the states at large momentum, where the band-mixing 
is negligible due to the marginal effective resonant cou- 
pling. Consequently these two optical conductivities are 
very close to the approximate formula Eq. (|C5[) in Ap- 
pendix [Cj However, at low Fermi energy, it is seen that 
although the dips due to the quasi-energy gaps also ap- 
pear in the optical conductivity from the projected dis- 
tribution, the peaks in the middle of the "steps" and 
the peaks from the states around zero momentum are 
absent. The absence of both effects is because the quasi- 
electron states with small k are occupied in the projected 
distribution. Thus the corresponding optical transitions 
are blocked. In addition, it is also shown that the van 
Hoff singularities do not affect on the optical conductiv- 
ity from the projected distribution, since the contribu- 
tion from the states around the momentum of the quasi- 
energy gaps is extremely small due to the strong band- 
mixing. 



IV. SUMMARY AND DISCUSSION 

In summary, we have performed a theoretical investiga- 
tion on the optical responses of graphene in the presence 
of intense circularly and linearly polarized THz helds via 
the Floquet theory. We examine the energy spectrum 
and DOS. It is found that gaps appear at small momen- 
tum in the quasi-energy spectrum, in consistence with 
the previous investigations ! 20 ' 21 ' 24 These gaps can be at- 
tributed to the ac Stark splittings induced by the single- 
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photon/multi-photon resonances. Nevertheless, in large 
momentum regime, where the energy spectrum has not 
been well investigated in previous works, we find that 
the quasi-energy gaps decrease dramatically with the in- 
crease of momentum and finally tend to be closed when 
the momentum is large enough. Consequently, taking 
account of the contribution from the states at large mo- 
mentum, the gaps in the DOS are effectively closed, in 
contrast to the prediction by Oka and Aokii^i 

We also investigate the optical conductivity from the 
mean-energy-determined distribution for different field 
strengths and Fermi energies. These results reveal the 
main features of the DFK effect in graphene. In the case 
with high Fermi energy, we discover that the optical con- 
ductivity presents a multi-step-like behavior around the 
optical frequency w; twice of the Fermi energy Ep, in 
contrast to the single-step-like behaviour in the field-free 
case. This effect is from the sideband-modulated optical 
transition, similar to the DFK effect in semiconductors. 
We also find that the optical conductivity varies mildly 
with the increase of u>i in each "step" , which is quite 
different from the DFK effect in semiconductors. This 
behaviour is due to the linear dispersion of graphene and 
the absence of the band-mixing at large momentum. 

The behaviour of the optical conductivity becomes 
more interesting in the case with low Fermi energy. We 
discover that dips appear at frequencies being the integer 
numbers of the applied THz field frequency f2, due to the 
quasi-energy gaps at small momentums. In the case with 
a circularly polarized THz field, it is found that peaks 
appear in the middle of the "steps" when the Fermi en- 
ergy is slightly lower than a local maximum of the mean 
energy. This kind of peaks become much less pronounced 
in the case with a linearly polarized field, owing to the 
anisotropic energy spectrum. Another interesting find- 
ing in the case with a circularly polarized field is that 
the contribution from the states around zero momentum 
can induce sharp peaks in the optical conductivity when 
the Fermi energy is lower than the mean energy at k = 0. 

Finally, we address the distribution function of the Flo- 
quet states. Our calculations are based on the ansatz 
that the distribution function of the Floquet states is 
determined by the mean energy, following the works in 
the literature! 28 ' 59 ' 60 The projected distribution function 
is also adopted in the literature. 26 It is noted that the 
multi-step-like behavior and the dips around frequencies 
IQ exist in optical conductivities from both distributions. 
This indicates that these two effects do not depend on the 
details of the distribution function and thus are expected 
to be observed in the optical absorption measurements 
subject to intense THz fields. Nevertheless, the peaks in 
the middle of the "steps" and the peaks from the states 
around zero momentum appear only in the optical con- 
ductivity from the mean-energy-determined distribution. 
By performing experimental investigation on these peaks 
in the optical conductivity, one can distinguish which dis- 
tribution function of the Floquet state is closest to the 
genuine one. One may also solve the kinetic equations 



with all the scattering explicitly included— to determine 
the distribution function. 
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Appendix A: Optical conductivity from extend 
Kubo formula 



We write Eq. (3) in Ref. [26| into the form 



A-7T 

OJl — 



X (fka - fk/3)6(£ ka - £k/3 + UJl) (Al) 

with f ka = J2 \{(ku\i/J ka )\ 2 n F (E ku ). Recall that E kv and 

V 

C, kv are the eigenvalues and eigenvectors of Ho. By using 
Eqs. flB) and (0, one obtains 



Re<r£?( W ,) 



4tt 



Ml *?2 



x - f£i)6(e km - e kt)2 -NC1 + Ul ) (A2) 

with = E^F(^)|(Ck,|^)l 2 - It is seen that the 

above equation is in the same form as Eq. (1311) . but the 
distribution f k ™ is quite different from f kr) = np(e kr i) 
used in the present paper. 



Appendix B: Approximate analytical solution of 
Schrodinger equation 

We first transform the effective Hamiltonian into the 
basis set formed by the eigenvectors of H$. Thus the 
Schrodinger equation can be written as 

1 d 
v F ka z + flrHz(t)] |*ta,(t)> = $^|*k,(*)). (Bl) 

Here \^ k ri(t)) — U^\^ kri (t)) , with the transformation ma- 
trix 



V2\ 1 



(B2) 



-ffrHz(#k,i) = U^HTHz(t)Uk can be divided into the in- 
traband H$[*(0 k ,t) and interband H^(8 k ,t) terms, 
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given by 



rrintra/-/) /\ 



printer//) .\ 
"THz (yk,*) 



/3f2<7 z (— cos cos #k sin fit 
+ sin #e sin 6>k cos fit) , (B3) 
/3f2<7 a (cos #e sin #k sin f2t 
+ sin6> E cos6>kCOsfH). (B4) 



Then we solve the Schrodinger equation without the in- 
terband term H^^*(9k,t) and obtain 



e- i " pfc V o) (0 k ,t)(l,O) T , (B5) 
e WFfct [u (o) (0 k ,t)]*(O,l) T , (B6) 



with 



u (°)(0 k ,i) = exp{i/3cos6> E cos#k[l - cos(fi*)] 

-i/3sm0 E sin0 k sin(f2t)}. (B7) 

It is evident that u(°\6\c,t) is a time-periodic function, 
thus the quasi-energy of l^knW) is exactly the same as 
the field- free energy rjvpk. This indicates that all quasi- 
energy gaps disappear without H^^(6k,t) ■ 

The next step is to write the solution of Eq. (|Blj) into 
the form 



i<m*)> = o,,i(*)i*Sw)+^. a (t)i*L-(*))- 

Substituting Eq. (|B8|) into Eq. (|B1[) . one obtains 



d t a, h2 (t) = -/3Qa^,i(t)^t/j(^k) 



3 -i(2i> F fc-Zfi)t 



(B8) 

(B9) 
(BIO) 



where 

2/; (#k) = (#k) (sin 6> E cos 6» k - « cos E sin 6> k ) 

+ w[|_\(^k)(sin6'E cos 0k + i cos#e sin#k), (Bll) 

J-o Jo 



a i2(3 cos #e cos 



Ji(2k(0k)|) 



-z{9v 



z(Q]i) = /3(sin6*E sin#k + i cos#e cos 0k)- 



(B12) 
(B13) 



In above derivation, we have applied the summation rule 
of the Bessel function 70 



with Z — x + e l9 y (x and y are real numbers). Since 
Eqs. (|B9[) and (IB10|) cannot be solved in analytical closed 
form, we solve these equations via the rotating-wave ap- 
proximation, which is widely used in the weak electro- 
magnetic field related problem.] 48 ' 49 Exploiting this ap- 
proximation, we neglect the rapidly varying terms with 
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FIG. 9: (Color online) Quasi-energy gap A m versus the cor- 
responding momentum from the exact calculation (dots) and 
the approximate formula Eq. (|B24[) (squares) under circu- 
larly polarized THz fields with /3 = 0.4 (solid curves) and 
2.3 (dashed curves). The triangles represent the results from 
the exact calculation without the intraband term flxSife- The 
indices m of the gaps at finite momentum from the exact cal- 
culation for p = 2.3 are labelled in the figure. The other 
gap indices can be read from the m-axis on the upper frame. 
Note that gaps only appear at the momentums where the dots 
(squares, triangles) indicate. The curves are only plotted as 
a guide for the eyes. 



/ 7^ m near the resonant point 2vpk ~ mQ and obtain 



1 



dta v ,i(t) = --pn arit2 (t)yl(e k ) 



Am* 



dta v ,i{t) = -PQa n)1 (t)y m (0 k )e 



—iSkmt 



(B15) 
(B16) 



with 5km — 2-upA; — mf2. Comparing the multi-photon 
resonance equations [Eqs. (|B15|) and (|B16[) ] with the 
well-known single-photon resonance equations in the two- 
level system [e.g., Eq. (83) in Ref. |48[, one can find that 
the main difference is the appearance of the coefficient 
Z/m($k)- 69 Thus \y m (0k)\ can be treated as a parameter 
to describe the magnitude of the effective coupling in- 
duced by the m-photon resonances. 

The solutions of Eqs. (|B15[) and (|B16|) are then easy 
to obtain: 



a n ,i = C+exp 



a v ,2 = C n exp 



it 
"2 



Y,j i0 JMJi-i(?)> (B14) Here 



(B17) 
(£18) 

(B19) 



A^^ k )=% ro ^ k )|0; 
C + and C'~ are the normalization coefficients satisfying 
|C+| 2 + |C-| 2 = 1, (B20) 



C+ 

Cn 



: S km + WSl m + [A' m (0 k )} 
pSly m (6w) 



(B21) 
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Thus one has 

|*k,(*)> - e-^ 5 lm+^'M? (c+e-i mnt u(°\e k , t), 
C-ei mnt [u^(9 k ,t)}*y. (B22) 
Evidently, the corresponding quasi-energy is 

^ = fV<*L + [AU0k)] 2 . (B23) 

Thus the quasi-energy gap at the resonant point 5k m = 
reads 

A ro (0 k )=/% ro (0 k )|tt (modfi). (B24) 

It is seen that the magnitude of the gap A m (0 k ) is de- 
termined by the effective coupling parameter |?/ m (# k )|. 

In Fig. [HI we plot the magnitude of the quasi-energy 
gaps A m against the corresponding momentums from 
the exact calculation (dots) and the approximate formula 
Eq. (IB24[) (squares) under circularly polarized THz fields 
with different field strengths. As shown in Figs. UK a) 
and (b), the momentums of the gaps from the calcula- 
tion markedly deviate from mil/ (2vp) in the strong field 
regime. Therefore, we define Ao as the gap at k = 
and A m (m 7^ 0) as the m-th gap from the one at k = 
for gaps from the exact calculation. The correspond- 
ing indices m are labelled in Fig. [9] In the case with 
low field strength (red solid curves in Fig. [9]), it is seen 
that the results from the approximate formula agree well 
with the ones from the exact computation. In the case 
with high field strength (blue dashed curves), the mag- 
nitude of the gaps from these two approaches are com- 
parable, but the corresponding momentums differ signif- 
icantly, owing to the rotating-wave approximation. In 
Fig. [9j we also plot the results from the exact calcula- 
tion without the intraband term i?™^ 3, (triangles). It 
is seen that except for the gap with m = 1, the gaps 
without -ff^H Z a are much smaller than the ones with the 
intraband term. In particular, the gaps with m being 
even numbers vanish within the error range of our com- 
putation (about 10~ 10 il). These results indicate that 
the intraband term plays a significant role in the 

formation of the quasi-energy gaps induced by the multi- 
photon resonances. 

Appendix C: Approximate analytical formula of 
optical conductivity 

Now we discuss the approximate analytical formula 
of the optical conductivity at high Fermi energy. It is 
known that the optical conductivity in this case is only 
from the contribution from the states with large momen- 
tum, where the effect from the interband term of H^Hz 
becomes negligible. Thus one can neglect the interband 
term and obtain the eigenvector 

\^ v (t)) = u k \yW(t)) (ci) 
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FIG. 10: (Color online) Time-averaged optical conductivities 
from the exact calculation (solid curves) and the the approxi- 
mate formula Eq. (|C5|) (dashed curves) are plotted as function 
of the optical frequency under a circularly polarized THz field 
with /3 — 2.3 for different Fermi energies. 

in which t/ k is given by Eq. dB2j) and \^{t)) is deter- 
mined by Eqs. (IB5|) and (|B6I) . Then one has 

EI07 = £ k»7 = r/v F k, (C2) 
= u^(e k )(e- ie -,l) T /V2, (C3) 
= [u W n (6 k )r(-e-^,lf/V2, (C4) 

where 

u (°)(0 k ) = oc. e*J n (\z(6 k )\) fc^f 

L \z{0u)\ . 

with z(0 k ) given by Eq. (|BT3|) . Substituting Eqs. flC2|) - 
([U5j) into Eq. ([3"Ij). one obtains 

ReojHwi) = 7- Y (! - N—)R N 0(u>i - 2E F - NO), 

(C5) 

where 

R N = d9 k Af{d k )J%(2\z{6 k )\), 
Jo 

with A x {0k) — sin# k and A y {9 k ) — cos6> k - 

In Fig. [TU1 we plot the optical conductivities from the 
exact calculation (solid curves) and the the approximate 
formula Eq. (|C5|) (dashed curves) as function of the fre- 
quency for different Fermi energies. It is seen that the dif- 
ference between two calculations becomes smaller when 
the Fermi energy increases. This indicates that the effect 
from the interband term of Hthz becomes negligible for 
high Fermi energy, in consistence with the discussion in 
the main text. 
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